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BONDED ORTHOTROPIC STRIPS WITH CRACKS* 

by 

F. Delale and F. Erdogan 
Lehigh University, Bethlehem, Pa. 

ABSTRACT 

In this paper the elastostatic problem for a nonhomogene- 
ous plane which consists of two sets of periodically arranged 
dissimilar orthotropic strips is considered. It is assumed 
that the plane contains a series of collinear cracks perpendi- 
cular to the interfaces and is loaded in tension away from and 
perpendicular to the cracks. First the problem of cracks fully 
imbedded into the homogeneous strips is considered. Then the 
singular behavior of the stresses for two special crack geome- 
tries is studied in some detail. The first is the case of a 
broken laminate in which the crack tips touch the interfaces. 

The second is the case of cracks crossing the interfaces. An 
interesting result found from the analysis of the latter which 
may have an important bearing on a possible delamination frac- 
ture Initiation at stress-free boundaries in bonded orthotropic 
materials is that for certain orthotropic material combinations 
the stress state at the point of intersection of a crack and an 
interface may be bounded whereas in isotropic materials at this 
point stresses are always singular. A number of numerical 
examples are worked out in order to separate the primary mater- 
ial parameters influencing the stress intensity factors and the 
powers of stress singularity, and to determine the trends regard- 
ing the influence of the secondary parameters. Finally, some 
numerical results are given for the stress intensity factors in 
certain basic crack geometries and for typical material combina- 
tions . 

This work was supported by NASA-Langley under the Grant NGR-39- 
007-011 and by the National Science Foundation under the Grant 
ENG77-19127. 


1 . INTRODUCTION 



In considering the failure a given structural component f 

if the corresponding material is homogeneous and isotropic in 
its strength and thermomechanical properties, the related frac- 
ture process is relatively well -understood and the techniques 
dealing with such problems are sufficiently wel 1 -developed. 

This is particularly true in the absence of large scale plastic 
deformations around the dominant flaw from which the fracture 
failure would develop. On the other hand in composites, parti- 
cularly in fiber-reinforced laminates, the situation is much 
more complicated not only because of the nonhomogeneity and 
anisotropy of the material which make it very difficult to ana- 
lyze the problem, but also because of the highly nonhomogeneous 
and nonisotropic distribution of the strength parameter making 
the development and the application of a proper fracture cri- 
terion also very difficult. In such materials it is quite 
possible that the concept of the progressive growth of a domi- 
nant crack with a well-defined leading edge is not an appro- 
priate model for the characterization of gross fracture beha- 
vior. Very often the damage zone developing around the dominant 
flaw is somewhat irregular and diffused and the fracture process 
is generally governed by a principle of "weakest link", the 
local fracture propagation being progressive or in discrete 
steps. Nonetheless, whatever the gross mechanism governing the 
process of fracture failure in the structure, one may nearly 
always assume that locally fracture initiation and propagation 
will take place along the leading edges of the existing flaws p 

where the conditions of the relevant fracture criterion are 
satisfied. Thus, in order to treat the local fracture pheno- 
menon in composite materials quantitatively, one may need the 
solution of the mechanics problem for flaws or cracks located 
at or near the phase boundaries or bimaterial interfaces. 

For composites which consist of bonded isotropic materials 
a wide variety of crack problems have been solved in which 
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either the asymptotic behavior of the stress state around the 
points of geometric singularity or the results for a specific 
crack geometry have been discussed (see, for example, [1] and 
[2] for review and references). Compared to the isotropic 
materials, the crack problems for homogeneous or nonhomogeneous 
anisotropic materials remain to be relatively unexplored. Most 
of the existing solutions refer to infinite planes [3-7]. The 
crack problem for an orthotropic strip is considered in [8] and 
that bonded to two orthotropic half planes is discussed in [9]. 
The details of the problem for a finite crack located in the 
neighborhood of or intersecting a bimaterial interface in 
bonded anisotropic materials do not seem to have been investi- 
gated, Even though the problem is rather complicated mostly 
because of the large number of independent constants entering 
the analysis, it may be managable under certain simplifying 
assumptions. The main assumptions made in this paper are (a) 
both materials are orthotropic, (b) the nonhomogeneous medium 
consists of two sets of periodically arranged dissimilar strips 
having different thicknesses, and (c) the cracks in the strips 
are collinear, perpendicular to the interfaces, and also periodi- 
cally arranged (Figure 1). Thus, one can take advantage of the 
symmetry of the medium and formulate the problem for two bonded 
strips only. The corresponding problem for isotropic layers 
or strips were considered in [10] and [11], In [12] the affect 
of the thickness and the elastic properties of the adhesive 
layer on the stress intensity factors in bonded dissimilar iso- 
tropic strips was considered. 

2. GENERAL FORMULATION OF THE PROBLEM 

Consider the plane problem for an orthotropic medium. 
Referring to, for example, [13] if u and v are the x and y com- 
ponents of the displacement vector, the equations of equili- 
brium may be expressed as follows: 
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for generalized plane stress, and 
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b-j-j/G^g * ^2 “ ^22^^12 > S3 ” 1 + *^12'^'"12 


(2) 


(3) 


for plane strain. Hera, v.^. , and G.^. , (i,o) = (1,2,3), 
are the engineering elastic constants, indexes (1,2,3) refer 
to the (x,y,z) directions, and the matrix (b^-j-) is given by 




(b,.g.) ■= B = A 


. A = (a,.j.) . (i,j) = (1.2,3) 


a..* = 1/E.. 1 a-. = -V../E.. = a 

The stress-displacement relations are 


, (i/j) 


„ = b M + b ^ 

^xx ^ll8x °123y 


„ - t. au . ' av 

* '^yy " “l 2 ax ° 22 ay 


(4) 


for plane strain, and 


lu = « 

ax Ej^j^ yy 

M + 11 „ fn 

By ax xy^ xy 


- n f BU . 8v \ 
^xy ~ ®xy^y ^ ax^ 


|v . V ^ ^ 

9y Eyy XX Ey^ 


(5) 


( 6 ) 


for generalized plane stress. 

Consider now the periodically arranged two sets of bonded 
orthotropic strips shown in Figure 1. In addition to the geo- 
metric symmetry indicated in the figure, it will be assumed 
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that the medium is loaded away from the crack region, parallel 
to the strips, and symmetrical with respect to the x-axis. 

Thus, the solution of the problem may be obtained by the stan- 
dard superposition technique, from the viewpoint of fracture the 
important component being the perturbation solution in which the 
crack surface tractions are the only external loads. One may 
note that because of symmetry it is sufficient to consider the 
problem for one quarter of each strip only. Let (x^ ,y) and 
(x 2 »y) be the local axes for the sets of strips 1 and 2 as 
shown in Figure 1. Let the displacements be expressed in terms 
of the following Fourier integrals: 

Vj(xj,y) = ijfj nij(a,Xj) sin ya da + ^ n^.(a,y) cos x^.a da , 


(7a, b) 

where j = 1, and j = 2 refer to the strips 1 and 2, respectively. 
Substituting from (7) into (1) one obtains a system of ordinary 
differential equations for the unknown functions f^ ,..,n^, which 
are coupled in pairs. Solving these equations we find 

^ S • r CtX - 

fj(a,x^.) = ^s_^Ajj^(a)e ^ ^ , mj(a,x-) * 



4* 

nj-(a,y) 


' • (j = 1.2) (8) 

In (8) Sjj, » (j=l ,2 , k=l,..,4) are the roots of the following 

characteristic equation; 
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The functions and B^.j^ 


1 (o=l)2 , k=l,..,4) are unknown. 


and the constants c^j^ and d^j^ are given by 


Cj-j " "C^.^ = (1 -B s )/B .,s^ 


'0 3 




*^3*2 " ■'^3*4 " 


^3l "" ■'*33 ^®jT^3‘l 

'*32 "'*3*4 “ ^^j2"^j|.^j5^^^03^j2^3*5 


( 10 ) 


The unknown functions Aj|^ and B^.j^ which appear in (8) are 
determined from the boundary and the continuity conditions of 
the problem. In addition to the assumed nature of symmetry in 
loading and geometry, it should be emphasized that in the per- 
turbation problem under consideration the only external loads 
are the local self-equilibrating crack surface tractions. 
Consequently, both components of the displacement vector would 


vanish for y-t-+w, and the x-component of the displacement, u^ 
{ 3 *= 1 » 2 ) would be zero along the axis of symmetry Xj = 0 , 

(3=1*2), Thus, the sixteen conditions which have to be used 


to determine the unknown functions A^^ and B.j^ (j=l,2; k=l,..,4) 


may be stated as follows: 


'3k 


UjVXjsVi-^-o , VjCxj,y)+o , (3*=1>2) for y-i-to , 


( 11 ) 


^ 1 (^ 1 *.^^) “ U 2 (“h 2 »y) , v*j(h-j,y) — 


0<y«» 

o-j j^^(h-( ,y) = '^2 xx^“^ 2’^^ ’ '^Ixy^*^! 


( 12 ) 


= cT 2 j^y{“h 2 *y) , 0<y<~ ,(13) 
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u^.(0,y) = 0 , = 0 > 0<y<» , (j=l,2) , (14) 


'^jxy(Xj’O) = 0 » , (j=l,2) , (15) 

” P-| V^-] ) > S 

V^(x^,0) = 0 , a<|x^i<h^ , (16a, b) 

^yyCxg.O) = P2(>^2^ » c<lx2l<d , 

V2(X2#0) = 0 , 0 <[x 2 !<c , d<]x 2 l<h 2 (17a, fa) 

In (9) It may arbitrarily be assumed that 

Re(Sj-j)>0 , Re (Sj. 2)^0 > (j=1»2) . (18) 

From (7), (8), (11) and (18) it therefore follows that 

Bj-^Ca) = 0 , 6^-2^“^ = 0 > (j=l,2) . (19) 


Ten of the remaining twelve unknown functions may be eliminated 
by using the homogeneous conditions (12-15) in (8), (7) and (6) 
The last two unknown functions are then determined from the 
mixed boundary conditions (ig) and (17), The problem may be 
reduced to a pair of integral equations by defining 

^Vj(xv,0) = <J>j(Xj) , 0<Ujl<hj. , (0=1,2) , (20) 

and by replacing the conditions (16) and (17) by (20). Thus 
all the unknown functions and may easily be expressed 
in terms of the new unknown functions ^ and We now observe 

that part of the mixed conditions, namely (16b) and (17b) is 
equivalent to 

■' a 

) “ 0 3 a'** [ x-j [ ‘'h*j 5 j 4*1 (x^ ) dx-j = 0 , ( 21 ) 

d 

‘^ 2 ^X 2 ) = 0 , 0<|x2l<c , d<lx2l<h2 ; |4>2(x2)dx2 = 0.(22) 

J c 

Through equations (8), (7), and (6), substituting the results 
into the conditions (16a) and (16b) we obtain two integral 
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equations to determine <|».j and ijjg. 

Because of the large number of elastic constants and un- 
known functions the process of deriving the integral equations 
is rather corapl icated and lengthy. However, the technique is 
straightforward and is quite similar to that followed in [10] 
and [11]. Therefore, the details of the derivations will not 
be given in this paper. As in [10], it can be shown that the 
integral equations are singular and may be expressed as follows 


I ClCt-x-j + t+x^ ^ kn(x-j»t) - (x^ ,-t)]^.j (t)dt 

^1 

t j" [k-j 2 ( »i* ) - k-j 2 ( » ” i*) ]<l2 ( ^ Pl(^l) > 


^1 3 


[kg^Cx^.t) - k2i(X2.-t)]4,,(t)dt + 





X 2 el .2 * [23a, fc 

where and Lg refer to the cracks on (y=0 , 0lx^<h-j) and 

(y=Q , 0<.Xg<h2) in the strips 1 and 2, respectively, and 

^1 " ^^lyy^l 4^ “'^l xy^lyx^ * ^2 ~ ^^2yy ^14^^^'^2xy^2yx^ 


C24] 

In deriving the integral equations one needs to define in a 
systematic fashion a large number of elastic constants and 
intermediate functions. Therefore, in order to conserve space 


the definitions leading to the expressions of the kernels 
k^-^. , (i , j=l ,2) , and the constants and y^^, and to the 


relationships between the functions A.^, 


and will also 


be omitted in this paper'' V . These definitions and the details 
of certain derivations may be found in [143 for the group of 
orthotropic materials which would give a characteristic equa- 
tion having only real roots Sj|^,(j = l,2; k=l,..,4) (defined 
henceforth as the orthotropic materials of type I), and in 
[15] for materials which would give a characteristic equation 
with only complex conjugate roots (defined, as the orthotropic 

f \ 

materials of type II)'- K 

The kernels k.,s which appear in (23) are of the following 

1 J 

form: 


k,^j(xTjt) — K^j(x*jtjci) 
■' 0 

Examining the behavior of K. . 


da , (i,j=l,2) . (25) 

for a-J-0 it can be shown that 



where c*- are known constants. 

■ J 

divergent kernels, by writing 


(26) 

Even though this may imply 



I, 


<|».(t)dt = 0 , (j=1.2) 


(28) 


it is s^en that the singularity at a = 0 may easily be removed. 


Also, by examining the behavior of the integrands (isj=l»2) 

for it can be shown that they decay exponentially provided 

the series of collinear cracks Lj (i=l,2) are fully imbedded in 
/★V * 

'Note that the constant T-ii is the same as the constants m,^ 
and r- 1 * defined in [8} (eqsl^l6 and 19) and the constants ‘ 
defined in (24) correspond to 4 p/(1+k) for the isotropic materials 

'■ 'in practice, since g.. in the characteristic equation (9) 
appears to be always a negative quantity, the third type of mater- 
ial giving four pure imaginary roots is not a realistic one. 
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the homogeneous strips (I.e., they do not touch or intersect 
the bimaterial interfaces). Thus, in solving the integral 
equations (23), may be treated as Fredholm kernels. In 
this problem since the kernels of the integral equations have 
only a Cauchy type singularity, the functions would have a 
square root singularity at the end points^'L-j and the equations 
may easily be solved by normalizing the intervals and by using 
the technique described, for example, in [16]. After solving 
the integral equations, the stress intensity factors may be 
obtained in terms of the functions tfi-j. For example, let Fig- 
ure 1 describe the crack geometry, i.e., let L-i = (0,a), 

(c,d); then, the stress intensity factors may be defined 
and obtained as follows [8]: 


k{a) = 


k(c) = 


k{d) = lim 
t-5-d 


Olyy(t.O) = -lim U-, 




>* 2 ’ 




4>2 ( t ) 


( i* ) 


(29a-c) 


3. CRACK TOUCHING THE INTERFACE 


Two limiting cases of the problem discussed in the pre- 
vious section are physically important and mathematically in- 
teresting. These are the cases of a broken laminate corre- 
sponding to a crack touching the interface (e.g., a=h-j, d<h 2 > 
Figure 1), and a crack intersecting the interface (e.g., a=h-| , 
d=hg, 0<c<h2> Figure 1). For example, referring to Figure 1, 
let a=h-| and d<h2. In this case it may be shown that as 
and for -hi<(xi , t)<hi , c<( jx 2 l > 1 1 1 )<d the integrands K-j 2 s ^ 2 -} , 
and K 22 in (25) decay exponentially. Therefore, the kernels 
^12* ^21* ^22 bounded in their respective closed do- 

mains. On the other hand for x-[-»-hi , t->-h'^ the exponential 
decay in K-|-,(xi,t, ) disappears, indicating that kn(xi,t) 
may contain terms which become singular as X] and t appr .^ch 
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tli6 end point hi simultaneously. These singular terms can be 
separated by studying the asymptotic behavior of the integrals 
given by (25) (see [16] for the technique and [10], [14], [17] 
and [18] for the application). To give an idea about the na- 
ture of these additional singular kernels let 

k^^(x^,t) = k-j^^(x.,,t) , 0 < (x^,t) < h., , 

(30) 

where k-jig represents the singular terms and is bounded 

in the related closed domain. Let the material be of type I 
with the real roots (see equation 9) 

s^1 = to-j>0 , S-J2 = “2^*^ ’ ^13 ~ “‘^l ’ ^14 "" '“2 ' 

(31) 


Then the asymptotic analysis of (25) would give 

(h-| -t)Bi5/tHi+h.|Ml 

(h-j-t)B-jg/m.| +u)2h-j 

[(h.|“t)e-[5/m-j-*-h-jU)2]'"f**'2''^l 

( h.| -t)S-'g/ui2^-hiiri-| 

( h*j "t) 6^1 g / 

C ( ^ 5 /^ 2 ^^ 'I ^2^ "(w22^'j) 

0 ^ (x^.t) < (32) 

where X 35 >..-> 3^88 constants and depend on the elas- 

tic properties of the materials only [14]. 
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Together with l/{t-xi), k-j-jg gives a generalized Cauchy 
kernel. Substituting from (30) into (23) and adopting the 
crack geometry shown in Figure 1 (with a=hi ) , the dominant 
part of (23) may be expressed as 

1 I- 1 

"J [-^73^+ 7 rk^-jg(x-[ ,t)](J)i (t)dt = (x-,) , -h-j<x^<h^ , 

-hi 

■] f d 1 

■? J ^2^ ^ ’^2 ‘"2- ’ c<Xg<d , (33a, b) 

where the bounded functions P-j and ?2 contain all the non- 
singular terms in (23). It is clear that the solution of 
(33b) is of the form 

^ 2 ^^) “ ’ c<t<d (34) 

giving the stress intensity factors as defined in (29). The 
singular behavior of the solution of (33a) may be studied by 
1 etting 


<i)l(t) = F^{t)/(hf-t^)^ , 0<Re(y)<l , -h^<t<h^ (35) 


and by using the function-theoretic method described in, for 
example, [16]. Thus, if we define the following sectionally 
holomorphic function 




1 fhl <f-](t) 

dt , (x^ = Re(z)) 


(36) 


by using (35) the asymptotic analysis of (36) gives 

1 Fl(h-j) ^ 




+ G ( z ) 3 

(2hi)^siniTY (z+h-j)^ (2hi)^simrY (z-h-j)^ ° 


lGo(z)l< 


I z+h^ 


Yo'^ReCy) 


(37a, b) 


where and are real constants. Now, substituting from 
(37) and (32) into (33a), we obtain the following characteris- 
tic equation to determine the unknown constant y: 


-Ecosiry + X 


“^1 ^ , ^1 ^ ^ ^2 ^ 

87^y^i-y 


M 1 *■ 

^2 




(38) 


where, again the material type I is assumed. It can be shown 
that for all material combinations (38) may have only one root 
in the strip 0<Re(y)<l and this root is always real. It can 
also be shown that as the orthotropic material constants tend 
to those for a pair of isotropic materials, the root obtained 
from (38) approaches the root of the corresponding isotropic 
characteristic equation given, for example, in [10] or [18], 

For this crack geometry the "stress intensity factor" may 
be defined in terms of the cleavage stress neigh- 

boring material which, from the fracture viewpoint, is the 
most important stress component. To calculate this we note 
that (23b) gives the expression for O2yy(x2s0) for -h 2 <.X 2 <.h 2 s 
that is in the uncracked as well as in the cracked portion 
of the strip. We also note that in the neighborhood of |x 2 |=h 2 
the singular behavior of cr 2 yy will be governed by the density 
function and the singular part of the kernel • As in K-j -j , 
it may be shown that for t-^-h-j , X 2 ->~-h 2 , the exponential decay 

in Kgi disappears, indicating that k 2 -j(x 2 ,t) may contain terms 
which become singular as X 2 and t go to the end point simul- 
taneously. If we again let 

^21 ^^2*^^ ~ ^21s^^2*^^ ^21f^^2*^^ * ( 3 ^) 

the singular part of the kernel may be separated and may be 
expressed as 
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• -I :■ 






(h|-'t)&|j/w^+a^hg 

{h-j-'t)p^g/u)T+a2hg 

^°^[(h^“t)p^g/u).j+(X2h2l^'(a25^2)^ 


^ ^ (^^T-t)g-jg/t02-^g^h2 

^*^^C(h-j-t)p^g/w2+ai hgl^-Ca-jXg)^ 

^°^C(hi-t)e-j5/w2+a2hg]"-(a2X2)^ 

Olhilh-j , Q^lXgllhg (40) 

where a-j and ag ai"e the positive roots Sg-j and Sgg of the 
characteristic equation (9) expressed for the strip 2 and the 
constants X are defined in [14], Thus, for the purpose of ana 
lyzing the singularity <^2yy expressed as 


*^2yy’^2'*^^ ~ ^2^ ^ ^21s^^2*^^*^l ^ ^20^^2^ (41) 

where pgg contains all the nonsingular terms. Upon substitut- 
ing from (35) and (40) into (41), the asymptotic analysis gives 


k(hn) 


(42) 


where *^q(^ 2^ remains bounded as Xg^-hg and the "stress intensity 
factor" k(h-]) is found to be 


k(h^) == (.t) 


„* - ^2 r:. , “Yr^^l ^1-Y 

^2 ■' Fin?r^^lOl“l ^^ 5 ^ ■^. ^ 102 “ 2 ^^ 3 ^ 5 ^ 


-Yf^2 W-Y . , 


104 “2 




( 43 ) 


4. CRACK CROSSING THE INTERFACE 


Consider now the case of a crack crossing the interface. 

In this problem the integral equations (23) are still valid 
with the two end points of the cuts L-j and Lg joining at the 
interface. For example, referring to Figure 1, let a = h-j, 
d = hgs and c>0. In this case at the end point x-j = h^ or 
Xg = hg all four kernels (x. ,t) will have singular terms. 
The singular parts k^-j^ and kg^^ coming from k^.j and kg^ are 
separated and are given by (32) and (40). Quite similar 
expressions may easily be obtained for k'j2s *^22s 
The dominant part of the system of singular integral equations 
may then be expressed as 


2 ( 

2, L. 

j=r J 


^iSl t"^ k^.j3(x^.,t)3(jij.(t)dt = Q^. (x^.) , 


= (-h^,h^) , Lg = (c,hg) , x.eL. , (i=l,2)(44) 

where in the analysis the symmetry condition of 02(Xg) = ■ 
-^>g("Xg) is used. If we now let 
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Cf ] 

FiCt) 

i * A i 

[t) - ^ 2 ^*’ 

r i' 

1 •] 




L . [. 




0<Re{e,S)<l j (45a, b) 

and define the following sectionally holomorphic functions 

h ^ c 

The asymptotic expressions for and Gg i'ay be obtained as [16] 
“ (2hi)esin^s (zH-h,)B ■ ’ 


62 ( 2 ) = 






(hr,"C)^simT6 (z-c)^ (h«-c)^simrp (z-h,)^ 


^5 n ( ^ ) > 


(47a, b) 


where G^.^(j = 1,2) has a behavior similar to that of 
which is given by (37b). Noting that outside their respective 
cuts G-j and Gg are holomorphic, substituting from (47) into 
(44), and following the procedure outlined, for example, in 
[16] (see, also [14] for details) we obtain 


F2(c)cotTr6 = 0 


Sf^-j(e)Fj.(hj) = 0 , (i=l ,2) , 


(48) 

(49) 


where the coefficients in the functions fj-(3) depend on the 

I u 

elastic constants of the two strips only and are given in [14]. 
Since Fg(c) and Fj(hj) , (a=1»R) are nonzero constants, (48) 
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gives the known result 6 = 1/2 and from (49) we obtain the 
following characteristic equation to determine the power of 
singularity 

aU) = = 0 , (i,0 = l,2) . 0<Re(g)<l . (50) 

It is also important to note that the end point values F^(h-|) 
and FgCh^) are not independent and are related by 

F2(h2) = -F^{h^)f^T(p)/f^2(p) (51) 

where e is the root of (50). An additional condition such as 
(51) is necessary to obtain a unique solution for the system 
of integral equations (23), since in this case there is only 
one singl e-val uedness condition which has to be satisfied by 
the displacement derivatives and namely 

[ ^(^g{t)dt ,+ 3>^(t)dt + (P2(t)dt = 0 . (52) 

■*0 -h*i "hg 

A systematic study of (50) indicates that for all material 
combinations the characteristic equations may have either no 
root or only a single real root in the strip 0<Re(p)<l. Also, 
p = 0 is always a root and there are no other roots with Re(p) 

= 0. In the foregoing analysis only the possibility of a power 
singularity is investigated. The results show that for certain 
material combinations (50) indeed has no root in 0<Re(p)<l, 
implying that for these materials at the intersection of the 
crack and the interface the stress state would be bounded. 
However, this analysis does not prove that in such cases there 
may not be a weaker, namely a logarithmic singularity. To 




investigate this question in (4.5) we let g = 0 and substitute 
the result into (46). We would then obtain the following 
asymptotic relations: 

F,(h,) Fa-h,) 

G.jU) = ■ - ■ M og(z~hi) ^p-^log{z+h^) + G^^(z) , 


SjCz) 


sinirs 


FyCh,) 

+ IT rlog(z-h„) 


+ Gg^Cz) , (53) 

where and Gg-j are bounded near and at the end points 
z = +hj and Gg^ has a behavior similar to (37b) in the neigh- 
borhood of z = c. Substituting now from (53) into the integral 
equations (44) we obtain 

F2(c)cotTr6 = 0 , (54) 

log(h^-x.)ig,jF^.(hj.) = R,(x,) , (1=1,2) (55) 

where and Rg are bounded functions and the constants g^.^ , 

(i,i=l,2) depend on the elastic constants only. Equation (54) 
again gives the known result 6 = 1/2. For (55) to be valid 
at X| = h- , (1=1,2) the coefficient of singular terms must 
vanish, or we must have 

2 g,;-* F.(h.) = 0 , (i=l,2) . (56) 

^ I J J U 

Since Fi(h.) , (j=l ,2) are nonzero constants, from (56) it 

J J 

follows that 
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To show (57) analytlcgny seems to be impossible. However, a 
systematic .numerical analysis indicates that for the material 
combinations having Q as the only acceptable power singu- 
larity (57) is indeed satisfied i identically. Furthermore, 
these studies also show that (56) always gives 


<1)1 (h-j ) (pf (hi ) 


^ -1 


(58) 


The result expressed by (58) meaning that in the composite 
medium the derivative of the crack surface displacement is con- 
tinuous at the interface Is, of course, the physically expected 
result. 

For the pair of materials in which (50) has a root in 
0<|i<l , at the poi nt (y=0, ^i“hi or X 2 =-h 2 ) the stress state 
will be singular. At this point, since the important stress 
■components are the normal and shear stresses on the interface, 
we may diTectly analyze the singular behavior of these stresses, 
to do this one has to go back to the original formulation of the 
problem and express these stresses in terms of the density 
functions «j>i and iig. Thus, after somewhat lengthy but straight- 
forward analysis we find [14,15] 

= i||^ ,(1=x.y) . 


by = (*"hi ,hi ) * ^2 ~ 
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studying the asymptotic behavior of the kernels h.^. it can be 
shown that as V'+O, s-t'+h-j tn<-'h|<;s<:h,j, and s-^-hg in .c<s<h2 simul- 
taneously the kernels become unbounded. By expressing 
h.^(y,s) = h^-j5(y,s) 

the singular parts h.. of these kernels can again be separated. 

I J S 

For example, for hj^^g(y,s) we obtain 

(hq+s)yn/2 (h-]-s)Yi/2 

(h^+s)^+(u-jy/&i5,)^ (h-j-sP+C^iJ^y/B^g)^ 

( h-j +s ) y^T’i 1 /2Yy2 ^ 1 1 ^^^1 2 

(hi+s.)^ + (ui2y/Pi5)^ (hy-sl^ + Cmgy/B-jg)^ 

, *^3^81 '*'^4^82 

^80 [{h,-s)BT5/to^.]=^ty^ 


V3^83‘**T4^84 th-j-s)g^ g/tog 


( 60 ) 


'^80 [ (h-| -3)8^ 5/ WgP+y- 

where the definition of the material constants y ^od X as well 
as the expressions for the remaining functions hj^2s> ^yls’ 


hy2s may be found in [ 14 ] and [ 15 ]. 


If the materials are such that the stress state at 
(y=0, x-[=hy) is singular, i.e,, 0<8<l , then one can again de- 
fine (fi. and S . , ( 1 = 1 , 2 ) as in ( 45 ) and ( 46 ) and obtain ( 47 ). 
Now observing that outside the cuts L-j and L2, specifically 
along the y-axis Gn and G^ are holomorphic, one can use ( 47 ) 







i-1 


i : 't 






to evaluate the singular terms in C59) (see, for example 
116^183). It can then be shown that 


^Ixx^h 


k..*: 


XX 




cf.(y) » (y>o) 


crixy(hi,y) = ^ + T^Cy) , y>o 


(61a, b) 


where the “stress intensity factors" may be expressed in terms 
of the density functions as follows: 


‘=xx “ Vxx ”J' • 


''xy - I'xy • 


(62a.b) 


The constants and y are known functions of the elastic 

A A Aj 


constants and may be found in [14] and [15]. 


For the material combinations in which and (j>g have 
no singularity at x-{ =" hi , Xg = -hg, (i.e., if 8 = 0 is the 
only acceptable root of (50)), since the kernels h^^ have 
singular pa>^ts of the form (60), from (59) it is not at all 
obvious that the stresses too would he bounded at the point 
(y=0, Xi=hi). This question can be examined by stustituting 
from (45) with B = 0, 6 = 1/2 into (59) and by going through 
a routine asymptotic analysis, which yields 
aixx(hi»y) = f^l i®xx ^ 

CTixy(hi»y) = D(y) (63a, b) 


i 


21 





where C(y) and D(y.) are hounded functions. It turns out that 
in all material combinations for which @ = 0 , the constant 
0 „\, is identically zero; therefore, the stresses are bounded. 

Considering the fact that in isotropic materials the 
stress state at the intersection of an interface and a crack 
is always singular (i.e., g> 0 ), from the viewpoint of delamina- 
tion or debonding fracture the practical importance of the pos- 
sibility of having bounded stresses at such locations in design- 
ing with certain orthotropic materials needs no elaboration. 

5. NUMERICAL SOLUTION 

In this paper the numerical results are obtained for 
several specific types of crack geometrieT In the first 
group of solutions it is assumed that the cracks are fully 
imbedded in homogeneous strips and [see Figure 1) 

a<h^ , c = 0 , d = b<h 2 • (64) 

The single crack, a = 0 , b 0 or a 0 b = 0 i s considered 
as a special case. In this problem the integral equations (23) 
are salved by using the Gauss-Chebyshev integration method [16] 
with L-j = {-a, a) , = (■'h,b) and under the single-val uedness 

conditions (28). The stress intensity factors are then obtained 
from (29a) and (29c) with d = b. 

In the second group of solutions it is assumed that a = h^ 
and 0 ^fa<h 2 - In this case the Gauss-Jacobi integration method 
is used to solve the integral equations. The details of the 
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numerical methoci may be found in 116,4 or 18J. After obtain- 
ing the density functions the stress intensity factors are 
calculated from f29c) (with d=^b) and (43). 

In the third group of solutions it is assumed that the 
crack crosses the interface, that is, a ~ h-[ s d = h 2 » o<c<hg 
(Figure 1). In this case for p>D, the integral equations (23) 
are solved by substituting from (45) and by using the Gauss- 
Jacobi integration method. Here the additional conditions are 
(51) and (52). After obtaining 4-] and stress intensity 

factors are determined from (29b) and (62) (see again [16], 
[14]1 or £18] for numerical procedure). 
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6. RESULTS 

The elastic properties of the materials used in the nu- 
merical examples are shown in Table 1. Materials 3, 4, and 6 
are basically isotropic and the remaining materials are ortho- 
tropic. For the materials 1 through 8 the roots of the charac- 
teristic equation (9) are real, meaning that they are of type 
I. Materials 9 and 10 are of type II for which (9) has complex 
conjugate roots. The numerical results given in this paper are 
all for the case of plane stress. Table H shows the material 
combinations used in the numerical analysis. The table also 
shows the powers of singularity y and p at the point of inter- 
section of the crack and the interface corresponding to a crack 
terminating at the interface (a=h-[, d<h 2 )> and that crossing the 
interface (a=h-], d^h^s oO), respectively (Figure 1). Unlike 
the isotropic materials, the characteristic equations (38) and 
(50) giving y and b in bonded orthotropic materials are quite 
complicated. They contain six independent material parameters 
and hence do not lend themselves to a relatively simple syste- 
matic parametric study. However, once the material combination 
is specified and g can be determined quite accurately. 

Even though it is very difficult to separate the material 
parameters which influence most of the values of B and y, and 
the stress intensity factors for the imbedded cracks, the cal- 
culations show that in this respect perhaps the most important 
single material parameter is the longitudinal stiffness ratio 
Eiyy/E 2 yy. lo Order to assess the effect of the remaining ma- 
terial constants a rather large number of calculations were 
done by fixing Eiyy and E2yy> systematically varying one at 
a time the remaining six constants, and by calculating y, Bj 
and the stress intensity factor k(a), the latter for imbedded 
cracks in material 1 only. The general trend is as follows: 

As E-|xx» ^ixy’ ^Ixy medium 1 containing the crack) 

are increased, y and k(a) increase, and as Egxx’ ®2xy* '’Zxy 
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Table 1 Elastic constants of the materials 
used in numerical calculations*. 


No . 

Exx/ 
109N/m2 
(1 O^psi ) 

Eyy/ 

1 O^N/m^ 
(1 O^psi ) 

Gxy/ 
109N/m2 
(lOSpsi ) 

vxy 

1(0) 

55,16 

(8.0) 

170.65 

(24.75) 

4.83 

(0.7) 

0.036 

2(0) 

134.45 

(19.5) 

31 .03 
(4.5) 

24.13 

(3-5) 

0.650 

3(1) 

154.77 

(22.447) 

155.83 

(22.6) 

59.68 

(8.655) 

0.300 

4(1) 

167.55 

(24.3) 

170.55 

(24.75) 

62.40 

(9.05) 

0.300 

5(0) 

10.07 

(1.46) 

31 .03 
1 (4.5) 

0.833 

(0.128) 

0.036 

6(1) 

30,34 

(4.4) 

31.03 

(4.5) 

10.83 

(1.57) 

0.400 

7(0) 

44.82 
1 (6.5) 

' 155.14 
! (22.5) 

4.83 

(0.7) 

0.020 

! 

8(0) 

34.48 
1 (5.0) 

1 6.895 

1 (I.P) 

3.45 
! (0.5) 

0.350 

9(0) 

21.37 

(3.1) 

i 66.88 
(9.7) 

17.93 

(2.6) 

1 

0.200 

10(0) 

1 17.24 

(2.5) 

17.24 

(2,5) 

6.895 

(1.0) 

1 ] 

i 0.760 

1 1 


(*) The materials are boron-epoxy and graphite-epoxy 
with various ply orientations. 





Table 2 The power of stress singularity y for a 
crack in medium 1 terminating at the in- 
terface and 0 for a crack crossing the 
interface. The properties of materials 
used in various combinations are given 
in Table 1 






Comb . 

Material s 

Power of Sing. j 

^lyy 

^2yy 


Med.l 

Med. 2 

Y 

3 

I 

1 

2 

0.55048 

0 

5 .50 

0.20 

II 

3 1 

9 

0.65699 

0.04248 

5.02 

2.48 

III 

4 

2 

0.66549 

0.04887 

5.50 

2.58 

IV * 

4 

6 

0.68914 

0.14547 

5.50 

5.75 

V 

4 

5 

0.80352 

0.05354 

5.50 

70.8 

VI 1 

7 

8 

0.74523 

0.05197 

22.5 

1 .40 

VII i 

2 

1 

0.42258 

0 

0.182 

' 5.00 

VIII 

2 


0.36911 

0.04248 

0.199 

0.403 

i 

IX 

9 


0.61554 

0.08520 

3.88 

i 2.6 

X 

10 

9 

0.43410 

0.08520 

0.268 

; 0.384 
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are increased, y and k(a) decrease. Among these variables 
the most significant factor influencing y and k(a) appears 
to be the ratio of shear moduli Gij^y/G2xy This may partly 
be observed also from Table 2 and Figure 2- The figure shows 
the stress intensity factor k(a) for imbedded cracks in ma- 
terial 1 as a function of the width ratio hg/h-j for a fixed 
relative crack length a/h-|=0.8 and for material combinations 

1, III, IV, and V given in Table 2. For these material pairs 
the stiffness ratio ^lyy/^Zyy constant whereas S-jxy/Ggxy 
is 0.2, 2,58, 5.75 and 70.8, respectively. It is seen that 
k(a} is consistently higher in material pairs having the 
greater Gixy/*=2xy i"atio. Figure 2 also shows that for h2-»-0, 
as expected, in all material combinations k(a) approaches the 
periodic collinear crack solution in an infinite plane which 
is the same for all homogeneous orthotropic as well as iso- 
tropic materials. 

A close examination of the results giving e, y, and kCa) 
indicates that generally one could accomplish a certain relax- 
ation in the stress singularity at the point of intersection 
of a crack and an interface in composites by introducing ortho- 
tropic materials. This may be seen, for example, by comparing 
the values for various material combinations given in Table 

2. In. fact for certain orthotropic material combinations it 
is even possible to have 0=0, i.e., no singularity, whereas 
in isotropic materials B is always positive, i.e., the stress 
state is always singular. The value of B has, of course, an 
important ^earing on the initiation of a possible delamination 
fracture from the stress-free boundaries . in bonded materials. 
Even though the result regarding the possibility of B=0 may 
appear to be somewhat paradoxial, considering the fact that 

in two isotropic wedges forming a half plane B is dependent 
on the wedge angles as well as the material constants and may 
be zero for certain ranges of wedge angles, it should not be 
completely unexpected. The possibility of reduction or com- 
plete elimination of singularity power B by varying the 
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secondary material constants seems to introduce an added 
flexibility in designing against the edge delamination in 
bonded materials. 

In solving the integral equations it is assumed that the 
composite medium is under a state of generalized plane stress 
and is subjected to external loads away from and perpendicular 
to the cracks. Thus the crack surface tractions in the per- 
turbation problem considered in. this paper are constant and are 
at the following ratio; 


Pi ^ = ^lyy 


Km 


(65) 


The stress intensity factors obtained for the imbedded cracks 
located in the first or second set of strips are given in Fig- 
ures 3-7. Comparison of the results given in Figures 3 and 4 
shows that for the same longitudinal stiffness ratio Eiyy/E2yy 
and the same material 2, k(a) calculated for an isotropic me- 
dium 1 is consistently greater than that calculated for an ortho- 
tropic material 1. This means that by introducing material 
orthotropy it is possible to obtain certain relaxation in the 
stress intensity factor. However, as seen from Figure 2, due 
to the effect of the secondary material parameters the opposite 
is also possible. In Figure 2 note that the combination IV 
refers to an isotropic-isotropic material pair - whereas III 
and V are isotropi c-ortho tropi c pairs giving stress intensity 
factors which are respectively lower and higher than that of 
IV, Corresponding results for the stress intensity factor k{b) 
for cracks imbedded in the second medium are given in Figures 
6 and 7. Materials in Figures 5 and 7 are of the type II and 
those in Figures 2, 3, 4 and 6 are of type I. It should be ob- 
served that as the thickness of the uncracked strips go to zero, 
the stress intensity factor in the cracked strips approach that 
of the periodic crack problem in the infinite homogeneous (iso- 
tropic or ortho tropic) medium. 
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Figure 8 shows a sample result for a composite medium in 
which both sets of strips contain cracks. Additional results 
for two as well as one set of cracks may be fou.nd in [14] and 
[15]. 

The stress intensity factors for the case of a broken 
laminate (i.e., for a=hi and c=d, or a=0 and d=h2» c=0) are 
given in Figures 9-12. Figures 9 and 10 show the results for 
material combinations I and II where all materials are of type 
I and Figures 11 and 12 give an example for the material com- 
bination IX where both materials are of type II. The figures 
show that in all cases as the width of the uncracked strip 
(i.e., the net ligament between the cracks) goes to zero, as 
expected, the stress intensity factors become unbounded. In 
these problems the stress intensity factor is defined by (42) 
and is calculated from (43). 

The results for a crack crossing the interface are given 
in Figures 13-19. In these problems the stress intensity fac- 
tor at the crack tip k(c) = k[j is defined by and calculated from 
(29b). For those material combinations in which b> 0 the stress 
intensity factors at the point of intersection of the crack and 
the interface kxx ^xy defined by (61) and are calculated 
from (62). For the material combinations II, IX and I used in 
these examples. Table 2 shows that power of stress singularity 
Y for a crack in material 1 touching the interface is greater 
than 1/2. Therefore, as the crack length 2 a approaches 2h-| or 
as c-^h2» the stress intensity factor kb at the crack tip calcuY 
Tated on the basis of 1/2 power becomes unbounded. Also, as the 
length of the net ligament 2c goes to zero kb again becomes un- 
bounded. These features of the solution may be observed from 
Figures 13, 16, and 19 giving the crack tip stress intensity 
factor as a function c/h 2 * Figures 13 and 16 show kb for ma- 
terial combinations II and IX in which 3>0. Figure 19 gives 
an example for the case in which p=Q. It may be noted that 
qualitatively the results for the two cases are quite similar. 
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The stress intensity factors and k^y for material 
combinations II and IX are given in Figures 14, 15, 17, and 
18. Note that in the limiting case of c^hg, that is for the 
case of the crack touching the interface, the power of the 
stress singularity at the interface would be y which is al- 
ways greater than s. Therefore, as expected and as seen from 
the figures, for c+h 2 the stress intensity factors calculated 
on the basis of singularity power p become unbounded. In 
these problems for the type of loading under consideration 
the normal component k^^^^ of the stress intensity factor seems 
to be negative. Since there is no crack surface interference, 
physically this means that normal stress along the interface 
near the crack surface is compressive, there is no inconsis- 
tency, and the singularity should be interpreted in the same 
way as in punch problems. 
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Figure 2 Stress Intensity factor k(a)=k^for cracks imbedded 
in strip 1 (0<a<h] , b-d=c) for material combinations 
I, III, IV and V, ^y~^lyy* ~^2yy» a/h]=0,8= 
constant ^ 
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Figure 12 Same as Figure 10, material combination IX 
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Figure H Normal component' kxx o'f the stress intensity factor 
at the intersection of the cracks and the in- 
terfaces in material combination II 
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